arXiv: 1506.0141 Ivl [cond-mat.supr-con] 3Jun2015 


Magnetic field dependence of the density of states in the multiband superconductor 

/5-Bi2Pd 


E. Herrera,^ 1. Guillamon,^’^ J.A. Galvis,^’^ A. Correa,^ A. Fente,^ R.F. Luccas,^ F. J. 

Mompean,^’^ M. Garcia-Hernandez,^’^ S. Vieira,^’^ J.P. Brison,^’^ and H. Suderow^’ 

^ Lahoratorio de Bajas Temperaturas, Departamento de Fisica de la Materia Condensada, 

Instituto Nicolas Cabrera and Condensed Matter Physics Center (IFIMAC), 

Universidad Autonoma de Madrid, E-28049 Madrid, Spain 
^Unidad Asociada de Bajas Temperaturas y Altos Campos Magneticos, UAM, CSIC, E-28049 Madrid, Spain 
^Departamento de ciencias naturales, Facultad de ingenieria, Universidad Central, Bogota, Colombia. 

"^Instituto de Ciencia de Materiales de Madrid, Consejo Superior de Investigaciones Cientificas (ICMM-CSIC), 

Sor Juana Ines de la Cruz 3, 28049 Madrid, Spain. 

^Univ. Crenoble Alpes, INAC-SPSMS, F-38000 Crenoble, France 
^CEA, INAC-SPSMS D-38000 Crenoble, France 
(Dated: June 5, 2015) 

We present very low temperature scanning tunneling microscopy (STM) experiments on single 
crystalline samples of the superconductor /3-Bi2Pd. We find a single fully isotropic superconducting 
gap. However, the magnetic field dependence of the intervortex density of states is higher than the 
one expected in a single gap superconductor, and the hexagonal vortex lattice is locked to the square 
atomic lattice. Such increase in the intervortex density of states and vortex lattice locking have been 
found in superconductors with multiple superconducting gaps and anisotropic Fermi surfaces. We 
compare the upper critical field Hc 2 {T) obtained in our sample with previous measurements and 
explain available data within multiband supercondutivity. We propose that /3-Bi2Pd is a single gap 
multiband superconductor. We anticipate that single gap multiband superconductivity can occur 
in other compounds with complex Fermi surfaces. 

PACS numbers: 74-25.-q,74.25Uv,74.55.+v 


INTRODUCTION 

Superconductivity is often found in binary metallic 
compounds with critical temperatures of the order of liq¬ 
uid helium temperature. Some are reviewed in Ref.[T] 
and give type II superconductors. Among them, MgB 2 
is peculiar, with a critical temperature Tc = 40K unsur¬ 
passed by related binary compounds. Such a high Tc re¬ 
sults from the combination of the strong electron-phonon 
coupling of the two-dimensional a bands and weak in¬ 
terband mixing with the three dimensional tt bands [2]- 
[8] . Each set of bands is derived from orthogonal orbital 
wavefunctions, leading to two well defined superconduct¬ 
ing gap features. Multigap and multiband superconduc¬ 
tivity have been conceptually linked together since the 
discovery of MgB 2 , suggesting that materials showing 
different Fermi surface sheets also have different super¬ 
conducting gaps in each sheet. 

Under magnetic fields, vortex core overlap is gov¬ 
erned by the Fermi surface velocity (see Refs. ISHn]). In 
MgB 2 , enhanced vortex core overlap has been observed 
and related to the Fermi surface properties of the sheet 
having a smaller sized superconducting gap[8l [T2IIT4] . 
This leads to a strong increase of the density of inter¬ 
vortex quasiparticle excitations measured by Scanning 
Tunneling Microscopy and Spectroscopy (STM) m US], 
and of the overall density of states measured by spe¬ 
cific heat or thermal conductivity, when applying a mag¬ 
netic field. The upper critical field shows a positive 


curvature, instead of the negative curvature expected 
within single gap s-wave BGS superconductivity pTl \Tb\ - 
|22]. The same occurs in many different compounds, 
including heavy fermions, borocarbides, and Fe based 
superconductors [TTJ l23ll29] . So far, the increased den¬ 
sity of states and positive curvature of Hc 2 close to Tc 
have been explained through multigap superconductiv¬ 
ity. However, the detailed mixed phase properties of a 
single gap s-wave superconductor and the possible influ¬ 
ence of multiband and of anisotropic Fermi surfaces in 
the mixed phase remain unclear. 


Here we show how superconducting features are modi¬ 
fied by a multiband Fermi surface. We measure /3—Bi 2 Pd 
(Tc = 5 K [m [30]) with a very low temperature STM. 
Our sample is moderately in the dirty limit (mean free 
path I smaller than superconducting coherence length 
£ < and the Fermi surface shows multiple sheets of 
mixed orbital character pT] [3T| (32 ] . The situation is op¬ 
posed to MgB 2 , with a and tt sheets that remain well 
separated even in presence of defectspJ)]. At 150 mK, we 
obtain atomic scale imaging and the hexagonal vortex 
lattice. At zero field we find single gap behavior follow¬ 
ing s-wave BGS theory but observe multiband properties 
in the mixed phase. The density of states in-between vor¬ 
tices increases more than expected for an isotropic sin¬ 
gle gap superconductor and the hexagonal vortex lattice 
locks to the square crystal lattice. 
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FIG. 1: Powder diffraction pattern of /3-Bi2Pd. Red symbols 
are the experimental points. The black line is the best fit 
to /3-Bi2Pd diffraction pattern [36]. Residuals are given by 
the blue line. The two series (upper and lower) of vertical 
green strikes represent, respectively, the position in 20 scale 
of the reflections from the yd-B^Pd (I4/mmm) and a-B^Pd 
(d2/ml) phases. Insets show a photograph of one /3-Bi2Pd 
crystal, and the temperature dependence of the resistivity. 


CRYSTAL GROWTH AND EXPERIMENTAL 
METHODS 


Single crystals of /3-Bi2Pd were grown using slight ex¬ 
cess of Bi [33l |34] . We grew our samples from high purity 
Bi (Alfa Aesar 99.99 %) and Pd (Alfa Aesar 99.95 %). Bi 
and Pd were introduced in quartz ampoules and sealed 
at 140 mbar of He gas. Then, ampoules were heated 
from room temperature to 900 °C in 3 h, maintained 24 
h at this temperature, slowly cooled down to 490 °C in 
96 h and finally cooled down to 395 °C in 200 h. This 
temperature is about 15 °C above the temperature for 
the formation of the a-Bi 2 Pd phase [35]. To avoid for¬ 
mation of the a phase, we quenched the crystals down 
to ambient temperature by immersion in cold water. We 
obtained large crystals of 5 mm x 5mm x 3mm. To 
characterize them, we made x-ray diffraction on crystals 
milled down to powder (Fig. 1, using x rays with wave¬ 
length 1.54 A). We find ;5-Bi2Pd (I4/mmm, see Rnf.[37]) 
with refined lattice parameters a = b = 3.36(8) A and 
c = 12.97(2) A and no trace of (a-Bi 2 Pd. We made in to¬ 
tal twelve growths, varying slightly the conditions for the 
quench, growth temperature and initial composition, and 
obtained always crystals with a resistivity vs temperature 
very similar to the one shown in Fig. 1. The tempera¬ 
ture dependence of the resistivity is shown in the inset 
of Fig. 1. The superconducting transition in our sample 
of /3-Bi2Pd crystals occurs at 5 K. Previous resistivity 
measurements in this material reported a slightly higher 
value of Tc (around 0.3 K larger) and a residual resistivity 
three time smaller than the one found here [17]. The spe¬ 
cific heat transition of our samples is sharp, of about 30 
mK width[38|, contrasting the transition width of about 



FIG. 2: a) Atomic scale imaging of the surface of /3-Bi2Pd 
at different length scales. Images are taken at a bias voltage 
of 10 mV and conductance 1/iS at 0.15 K. Lattice structure, 
highlighting the cleaving plane with distance between two suc¬ 
cessive planes d = 6.6A is also shown. The inset in the top left 
image shows a cut through a line showing a jump whose size 
corresponds to the distance between cleaving planes. Bottom 
right panel is a Fourier transform of the atomic size images. 
Arrows give crystalline axis, b) Temperature dependence of 
the experimental tunneling conductance (black dots). Red 
lines are fits to the s-wave BGS expression at each temper¬ 
ature, leaving A as the only free parameter. The values of 
A obtained are plotted in the bottom panel, together with 
the temperature dependence obtained from BGS theory (red 
line). 


300 mK reported in Ref. m- We measured also the upper 
critical field Hc 2 {T) using resistivity and susceptibility as 
a function of temperature or magnetic field for the field 
applied parallel to the c-axis. The results coincide with 
the positions where we also observed vanishing supercon¬ 
ducting features in STM tunneling conductance. 


To make the STM measurements, we use a home built 
set-up installed in a dilution refrigerator with an energy 
resolution in the tunneling spectroscopy of 0.15 K. Con¬ 
struction is similar to Ref. [39]. We also took a few tun¬ 
neling data at fields parallel to the surface using a three 
axis coil system described in Ref. [40]. We use an An tip 
cleaned by repeated indentation on an An sample as de¬ 
scribed in Ref. mi- We make the STM measurements in 
a sample roughly one mm thick, which was cleaved using 
a scotch tape at ambient conditions after glueing it using 
silver epoxy to the sample holder. Usually, bias voltage 
is kept at 10 mV or below, and the tunneling conduc¬ 
tance is of a few tenths of fiS. Topography and vortex 
lattice images are independent of tunneling parameters. 
To obtain vortex lattice images we cut the feedback loop 
at each point and make full I-V curves, as in previous 
wor kiia. No filtering or image treatment is applied to 
the topography and conductance maps shown here. 
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FIG. 3: The vortex lattice as a function of the magnetic held applied parallel to the c axis direction of /3-Bi2Pd. The vortex 
lattice has a hexagonal symmetry for all values of the magnetic held. The images were taken at 150 mK and the magnetic 
held was increased up to Hc 2 =0.6 T. Constrast in the zero bias conductance is shown with color scales. The orientation of the 
crystalline axis is shown as green and blue arrows in each panel. 


RESULTS 

The tetragonal structure with Bi-Pd blocks (Fig. 2) 
suggests that it is easy to obtain clean and atomically 
hat surfaces by cleavage or exfoliation. In Fig. 2a we 
show atomic resolution topography images taken on the 
surface of the sample at 0.15 K. We hnd indeed a square 
atomic lattice. The Fourier transform of the topography 
images gives a lattice parameter of a = b = 3.3 A, coin¬ 
ciding with crystal structure values. Surfaces are atomi¬ 
cally hat over hundreds of nm. Small steps are sometimes 
viewed in the images. The step in Fig. 2a (top left panel 
and inset) is of 6.5 A height, which corresponds to the dis¬ 
tance between adjacent Bi-Pd groups (black arrow in top 
right panel of Fig. 2a). Bi-Pd groups are more strongly 
coupled than the Bi-Bi sheets, because the directional 
Bi-Bi bonds are weaker than the Bi-Pd bonds[32]. Thus, 
we conclude that the surfaces in Fig. 2a are made out of 
the square Bi lattice. 

Fig. 2b shows the tunneling conductance vs bias volt¬ 
age as a function of temperature. At 0.15 K, we hnd 
clear superconducting quasiparticles peaks and no con¬ 
ductance at zero and low bias. We can ht our data using 
single gap BCS theory and A = 0.76 meV. Temperature 
dependence of the superconducting features shows that 
superconductivity disappears at about 5 K. 

When we apply a magnetic held parallel to the c axis, 
we observe an ordered hexagonal vortex lattice on large 
atomically hat regions. In Fig. 3 we show vortex lattice 
images from 5 mT to 500 mT. A hexagonal Abrikosov 
lattice is observed in all images. The intervortex distance 
is modihed as expected for the hexagonal vortex lattice 

= 1.075a( with 00 being the hux quantum). 


We hnd that the orientation of the vortex lattice is 
determined by the underlying crystalline lattice (arrows 
in Fig. 3). One of the three main vortex lattice directions 
is always parallel to one of the two crystalline axes. This 
gives two equivalent orientations for the hexagonal vortex 
lattice at any magnetic held. 

At a hxed magnetic held, we hnd the two different 
vortex lattice orientations when changing the scanning 
windows (Fig. 4), obtained by moving in-situ the sample 
holder using the method described in Ref. [39]. Thus, the 
orientation of the hexagonal vortex lattice forms domains 
oriented along one crystal axis. The size of these domains 
is considerably larger than the scanning window (2 jam 
X 2 /im). 

The spatial dependence of the superconducting den¬ 
sity of states in and around vortices shows that the 
gap fully closes inside the vortex core (Fig. 5a). We do 
not observe signatures of Caroli-de Gennes-Matricon An¬ 
dreev core states [43ti45] . From the residual resistivity 
of our samples (p = 18 iiVLcm just above Tc) we esti¬ 
mate the mean free path using Drude formula and hnd 
^ = 15.3 nm. On the other hand, the in-plane coher¬ 
ence length from the upper critical held (discussed below, 
Hc 2 (F = f)K) = 2 ff 2 ) yields ^ = 23 nm. Thus, ^ ^ and 

Caroli-de Gennes-Matricon states are smeared by defect 
scattering [46], 

In Fig. 5b we show the spatial dependence of the zero 
bias conductance cro(r) along a path crossing several vor¬ 
tex cores at 0.3 T in 0—Bi 2 Pd (blue points), and the 
expectation assuming vortex core overlap far from Hc 2 - 
The latter is calculated by summing over relevant neigh¬ 
bors (r0 using cro(r) = 1 — tan/i((r^ — r)/^) (dark 

dashed line) with ^ = 23 nm. This is an approxima¬ 
tion widely used in literature, roughly conhrmed by mi- 



4 



FIG. 4: Vortex lattice at three different magnetic fields. In 
each magnetic field we show two vortex images obtained at 
different scanning windows. Note the change in the orienta¬ 
tion with respect to the atomic lattice, as marked by green 
and blue arrows. Contrast (zero bias conductance) is shown 
with the color scales. 

croscopic calculations which show that the intervortex 
density of states for single band superconductors in the 
dirty limit is indeed practically negligible for fields below 
about half Hc 2 |47H49] . Our data show that the tunneling 
conductance in-between vortices is more affected by the 
magnetic field than expectations for a single gap super¬ 
conductor. The magnetic field increase of the intervortex 
density of states is pronounced (open blue points in Fig. 
5c and Fig. 6 ). It is smaller than the increase found in 
the superconductors 2 H-NbSe 2 , 2 H-NbS 2 and MgB 2 but 
also above the increase expected for a single gap s-wave 
superconductor (dark dashed line in Fig. 5c). When we 
apply the magnetic field along the basal plane, we find a 
smaller sized intervortex density of states (Fig. 5c). 

We have measured Hc 2 (T) using resistivity and suscep¬ 
tibility (Fig. 7), and compared results with available data 
in a sample with a larger mean free path [17]. A smaller 
mean free path leads to a shorter coherence length and 
hence to an increased Hc 2 (T), also in multigap supercon¬ 
ductors, see for instance MgEofTOl [501453] . Here, how¬ 


ever, we do not observe such an increase in Hc 2 and in¬ 
stead find the same result as in previous measurements 
on samples with a larger mean free path. 

To explain this result, we have calculated Hc 2 {T) us¬ 
ing a multiband approach and the method described in 
Refs. p!5l[Tb][T8142Q] (see also Annex for details). We sim¬ 
plify the Fermi surface by using a two band description, 
with differing Fermi velocities VF,i^ electron-phonon cou¬ 
pling Xij and electron scattering parameters Tij (with 
subindices i, j for each band). We assume that the intra¬ 
band and interband electron-phonon coupling constants 
Xij are connected together as An = A 21 and A 22 = A 12 . 
The superconducting gap has then the same value in 
both bands. This reproduces our zero field tunneling 
spectroscopy result (Fig. 2 ). Specific heat measurements 
made very recently in our samples are also compatible 
with a single superconducting gap [38]. For the electronic 
specific heat coefficient we use 7=13 mJ/K‘^mol and de¬ 
rive the partial densities of states for each band from 
the anisotropy in the Fermi velocity. We also correlate 
the values of the Tij with the experimental values of the 
resistivities. To this end, we assume that the conduc¬ 
tivities of each band are equal pi = p 2 and we deduce 
the relaxation rate for each part of the Fermi surface 

2 

using Ti = . Where Vmoi is the molar volume, 

7 the Sommerfeld coefficient of the specific heat and p 
the resistivity. Moreover, A 12 /A 21 is given by the ra¬ 
tio of the density of states in both bands. Thus, taking 
7 = 13 we deduce 7 ^ from the choice of the 

Xi^j. We also assumed that the relaxation rate is isotropic 
and that each is split into intraband and inter band 
scattering, parametrized by an additional coefficient, a, 
such that Til = QfTi, ri 2 = (1 — a)Ti. The introduction 
of interband scattering ri 2 7 ^ 0 is needed to decrease the 
sensitivity of Hc 2 (T) to the mean free path ^ and obtain 
similar values of Hc 2 {T) when i is decreased. In table I, 
we give the list of parameters values used to fit the two 
sets of data. 


We neglect interband defect scattering for the data of 
Ref. [T7| . This provides an excellent fit to the Hc 2 {T) 
data of Ref. [17] (Fig. 7). If we simply decrease the mean 
free path to try to fit Hc 2 (T) in our sample, we find an 
increase of the upper critical field. Only when allowing 
for interband scattering, with Tij 7 ^ 0 for i, j 7 ^ 0 , we find 
that Hc 2 (T) does not increase when the mean free path is 
decreased. An intraband decrease of the mean free path 
invariably leads to a decreased coherence length and an 
increased upper critical field, but interband mixing can 
lead to similar values of the upper critical field. Intro¬ 
ducing interband scattering, we obtain an excellent fit of 
our data, with, in particular, an upper critical field that 
does not increase with decreased mean free path(Fig. 7). 
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FIG. 5: a) The tunneling conductance vs bias voltage when entering a vortex core at 0.3 T along the path shown by the white 
arrow in the inset (one curve each 6 nm approximately). In b) we show the normalized zero bias conductance cro as a function 
of the position at 150 mK and at i7/i/c2 = 0.5 (left panel, data shown as points, continuous blue line is a guide to the eye and 
black dashed line is calculated as discussed in the text), c) we show the magnetic held dependence of the zero bias tunneling 
conductance cTmid exactly in-between vortices (blue points, continuous blue line is a guide to the eye) when the held is applied 
along the c-axis (open blue points) and when the held is applied along the plane (hlled blue point). Data from 2H-NbSe2, 
2H-NbS2 and MgB 2 (points, lines are guides to the eye) are Refs. [T^ fT^ W2[ HS] . 


Parameters for the calculation of Hc 2 (T) 


Sample 

VFl 

VF2 7 


mol) p (pQ, cm) 

Tc (K) 

a P 

This work 

0.09 

0.4 

13 

18 

5.07 

0.2 0.8 

Ref. Imai [17] 

0.09 

0.4 

13 

5 

5.35 

1 1 


TABLE I: Parameters used to calculate the temperature dependence of the upper critical held along the c-axis in samples used 
in this work and in samples of Ref. [17] (Fermi velocities are given in units of 10^ m/s). 


DISCUSSION 


The inhuence on the vortex lattice of a square crystal 
symmetry has been discussed in detail before. At low 
magnetic helds, neutron scattering studies of the vortex 
lattice in tetragonal superconductors such as TmNi 2 B 2 C 
and CeCoIus show two vortex lattice domains [131 154H56] . 
In other tetragonal nickel borocarbides, and in VsSi, the 
square symmetry produces a transition between hexago¬ 
nal and square vortex lattices when increasing the mag¬ 
netic held above helds of the order of a Tesla. The 


current distribution around vortices is sensitive to non¬ 
local effects, which introduce a radial dependence in the 
vortex-vortex repulsion due to the shape of the Fermi 
surface [54l (551 [57f[59] . In /3—Bi 2 Pd, the square electronic 
symmetry does not give transitions in the vortex lattice 
symmetry. But the vortex lattice is still locked to the 
crystal lattice. Thus, non-local effects remain, showing 
that the Fermi surface features are playing a prominent 
role in the orientation of the vortex lattice. 

The lack of variation of Hc 2 with mean free path points 
towards the inhuence of multiband Fermi surface on the 
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FIG. 7: a) Hc 2 (F) along the c-axis in our sample of /^-BbPd 
with a residual resistivity of 18 /iQcm (open blue points) and 
Hc 2 (F) in sample of Ref.|17j with a residual resistivity of 5 
/iQcm (open red points). Lines are fits to each Hc 2 (F) using 
the parameters explained in the text. Note that the upper 
critical field remains nearly the same, although the mean free 
path varies by a factor of four. Magnetic field is applied along 
the c-axis. In the inset we show the transition in resistivity 
and susceptibility in our sample at a magnetic field of 0.1 T. 
The transitions remain sharp and well defined throughout the 
phase diagram. 


FIG. 6: Magnetic field dependence of the full normalized tun¬ 
neling conductance cTmidiV) obtained exactly at the center 
between vortex cores. Magnetic field is applied along the c- 
axis. Gurves are shifted vertically for clarity. 


mixed state properties. In MgB 2 , two-dimensional sheets 
with strong electron-phonon coupling are derived from 
the a electrons of B orbitals, whereas the three dimen¬ 
sional sheets are derived from tt band orbitals. Inter¬ 
band scattering is particularly small in MgB 2 . Previ¬ 
ous theoretical work remarked that a decrease in Hc 2 
might be observed by producing strong enough inter¬ 
band scattering [TOlHO]. The experiments show, however, 
a strong decrease of Tc with interband scattering, con¬ 
comitant with a more isotropic superconducting gap and 
increased Hc 2 [ 5 TH 53 I [60] . 

Bi 2 Pd has a Fermi surface with two nearly cylindri¬ 
cal sheets and two 3D sheets, derived from Pd 4d and Bi 
6p states [3T1|32]. States from both Pd and Bi contribute 
to the density of states at the Fermi level, and in partic¬ 
ular anisotropic Adxy+yz orbitals. The strong anisotropy 
of these orbitals favors interband scattering. With a sin¬ 
gle superconducting gap, the zero field Tc is not strongly 
affected by interband scattering, yet Hc 2 shows the pe¬ 
culiar behavior discussed here. 

We now compare our results results on the deviation 
in the magnetic field dependence of intervortex density 
of states with respect to expectation for single gap super¬ 



FS parameters from Hc 2 {T) 

A (meV) 


vfi (10® m/s) vf2 (10® m/s) 

V 

<1 

A 

<1 

MgB2 

0.29 

0.9 

7.1 2.2 

2H-NbS2 

0.155 

3.1 

0.97 0.53 

2H-NbSe2 

0.055 

1 

1.2 0.75 

/3-Bi2Pd 

0.09 

0.4 

0.75 0.75 


TABLE II: Parameters used to account for Hc 2 {T) of MgB 2 , 
2H-NbS2, 2H-NbSe2 and Bi 2 Pd. The gap values A> and 
A< are obtained from [TJl [131 SH SSI [60] . 


conductors in /d-Bi 2 Pd with known features of MgB 2 , 2H- 
NbSe 2 and 2H-NbS2 (table II) [H [El [IHHZQI 021 SSI |60] . 
MgB 2 is the compound where the two gaps are more sep¬ 
arated in energy, being the larger gap a factor of three 
higher than the smaller gap. Moreover, the gap distri¬ 
bution is narrow around these two values providing two 
neat features in the energy dependence of the supercon¬ 
ducting density of states [60]. The latter is also found 
in 2 H-NbS 2 although the ratio between the two gaps is 
somewhat smaller, around two [42]. The same gap ratio 
is found in 2 H-NbSe 2 , although in this material there 
is a sizeable in-plane gap anisotropy with a wide distri¬ 
bution of gap values |43l [CTH63] . This makes the vortex 
core overlap in 2H-NbSe2 (Fig. 5) smaller than in 2H- 
NbS 2 . On the other hand, the anisotropy in the Fermi 
velocity is stronger in 2H-NbS2 and 2 H-NbSe 2 than in 
MgB 2 m [13 [nmni sa sa . Yet, the vortex core overlap 
is the highest in the latter. Thus, the strongest increase 
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in vortex core overlap is produced by multigap supercon¬ 
ductivity. 

In conclusion, we have obtained atomically flat Bi 
surfaces in Bi 2 Pd where we observe an isotropic su¬ 
perconducting gap and a hexagonal vortex lattice. By 
discussing tunneling spectroscopy, vortex lattice and 
Hc2 (7"), we have shown that yd—Bi 2 Pd is a multiband su¬ 
perconductor with a single superconducting gap. Inter¬ 
band scattering precludes the usual increase of Hc 2 with 
a decreased mean free path. The hexagonal vortex lattice 
orientation locks to the crystalline lattice. We conclude 
that the mixed phase of superconductors is strongly mod¬ 
ified in multiband Fermi surface materials, even when the 
zero field superconducting density of states is not. 
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ANNEX I. CALCULATION OF THE UPPER 
CRITICAL FIELD. 

Several papers provide methods to calculate the upper 
critical field Hc 2 (T) in multiband superconductors pdH[TSl 
[T6l[l8H20l[50l[g. Here we follow Ref. m , a microscopic 
calculation of Hc 2 {T) taking into account defect scatter¬ 
ing. The upper critical field Hc 2 {T) is found by calcu¬ 
lating the set of Pi with largest values that solve the 
following equations: 

'[Ui{n) = uJn^TrT - n) + 

j,m 

Smn {nj / 27rT)) sgn('a;^) (1) 


Ai{n) = ttT - n) - - \^m\) + 

j,m 

^mn {^i,j /27rr)]A:’j(m)Aj(m) (2) 


Xi{n) = {2/\fPi) / dqe ‘^^tan ^{q^J~p^/{\wi{n)\ + 

Jo 

YMB-f^c2Sgn(w„))(3) 


A = (4) 

KM)= ! ( 5 ) 

ujn are the Matsubara frequencies, uja‘fjF{uj) is the 
electron-phonon coupling and the Cooper pair wave- 
function, Tij the relaxation rate, Xij the electron phonon 
coupling constant and vp^i the Fermi velocity at a plane 
perpendicular to the magnetic field in each band i. The 
Fermi velocities discussed here are unrenormalized, i.e. 
as obtained without the pairing interactions. Their val¬ 
ues eventually found in experiments (as quantum oscil¬ 
lation or photoemission) need to be renormalized by the 
electron-phonon interaction Xij. 

We see that the relevant parameters to describe mix¬ 
ture between different bands are the off-diagonal terms 
of matrices Xij and Tij. The superconducting order pa¬ 
rameter is found by an equation of the same form of the 
BCS self-consistency gap equation (2). The first term of 
the equation between square brackets [ ] accounts for in¬ 
traband and interband scattering due to electron-phonon 
interaction and defect scattering. 

The bare Fermi velocity vp^i enters into the equation 
through the term Xj , which also depends on the electron- 
phonon coupling and the interband and intraband scat¬ 
tering parameters Xij and through equation (1). 
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